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ABSTRACT 

This  thesis  is  concerned  with  proving  that  the  eigen- 
values of  a  specific  unsymmetric  matrix  are  real  and  positive, 
without  actually  computing  them.   The  method  of  finite 
differences  is  applied  to  the  vibration  analysis  of  a 
cantilever  beam  and  leads  to  an  unsymmetric  stiffness  matrix 
in  the  eigenvalue  problem  formulation.   The  technique 
employed  in  the  proof  is  based  on  a  perturbation  theory 
given  by  Wilkinson  for  real  symmetric  matrices.   Application 
of  the  theory  is  made  to  the  cantilever  beam  eigenvalue 
problem o   The  results  verify  that  the  eigenvalues  of  this 
and  other  unsymmetric  matrices  can  be  proven  real  and 
positive  without  their  actual  values  being  calculated. 
A  multiple-variable  formulation  of  the  cantilever 
beam  vibration  analysis  is  also  examined  to  illustrate  that 
it  implicitly  contains  all  of  the  properties  of  a  symmetric 
matrix. 
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I.   INTRODUCTION 

The  vibration  analysis  of  many  structures  can  be 
accomplished  by  the  application  of  finite  differences  to 
the  partial  differential  equations  governing  the  free 
undamped  vibration  of  the  structure.   Assuming  the  motion 
of  the  structure  to  be  harmonic  _,  with  natural  frequency  CO  s 
leads  to  the  matrix  equation  for  the  eigenvalue  problem 

2 
Kc56=CO  M0  (1-1) 

where  K  and  M  are  square  matrices  and  are  referred  to  herein 
as  the  stiffness  matrix  and  the  mass  matrix  respectively. 

symmetric,  and  K  or  M  is  positive  definite ,  the  eigenvalues, 

2 
CO    5    are  real;  when  both  K  and  M  are  positive  definite , 

CO      is  always  positive   Ref.  1   .   When  K  and  M  are  developed 

using  energy  techniques,  they  are  always  symmetric;  however, 

when  they  are  developed  for  the  same  problem  by  applying 

finite  differences  to  the  governing  partial  differential 

equations  of  motion  and  the  boundary  conditions,  they 

frequently  emerge  numerically  as  unsymmetric,  calling  into 

2 
question  whether  numerical  results  for  oj   will  in  fact  be 

real. 

In  this  thesis  an  unsymmetric  stiffness  matrix  is 

developed  by  finite  differences  for  the  cantilever  beam. 

A  perturbation  theory  for  real  symmetric  matrices  given 


by  Wilkinson   Ref.  2   is  used  to  establish  that  all  of 
the  eigenvalues  of  the  unsymmetric  matrix  formulation  are 
real  and  positive;  hence  the  results  obtained  for  the 
natural  frequencies  are  real^,  representing  undamped  harmonic 
motion  as  assumed 0 
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II.   DERIVATION  OF  THE  MATRIX  EQUATIONS 

The  structure  selected  for  analysis  is  the  cantilever 
beam  of  length  L  shown  in  Figure  1. 
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BEAM     CONFIGURATION 
FIGURE  I 


-►X 


The  equation  governing  the  free  undamped  vibration  of  the 


beam  can  be  given  in  the  form 


r$ 


H- 


n2 


(2-1) 


in  which  t  is  time,  and 

$=  Y/L 
x  =  X/L 
/LI  =  p/(Er2) 
where  Y  is  the  deflection  of  the  beam,  X  is  the  axial 
coordinate,  a  is  the  mass  density,  E  is  the  elastic  modulus, 
and  r  is  the  radius  of  gyration  of  the  cross-section. 
Assuming  harmonic  motion, 

<f>  =c±>  sin  art  (2-2) 
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Substituting  Equation  2-2  into  Equation  2-1  leads  to 

4%  -M>  (2-3) 

d  x 

2 
in  which    \=UCO  . 

The  boundary  conditions  on  the  beam  are 

^  _  deb  A 

<P  ~^~~   °      at  x  =  o  (2-4a) 

and 

d2^  =  d3cfr_  0 

dx2    dx3  at  x  =  i  (2-4b) 

Let  the  beam  be  divided  into  ten  equal  intervals  of 
length  hL  as  shown  in  Figure  1.   The  nodes  j=-l5  0=11.?  and 
j=12  are  fictitious  nodes  located  off  of  the  beam.   Let 
the  derivatives  in  Equation  2-3  and  Equations  2-4  be 
approximated  bv 


*&\     -1      (  -0^  +  0_..J  /  (2h)  (2-5a) 

J 


/ 


( 

(^f\  -      (  0,  ,  -  2«,  +  0.   )  /  h2     (2-5b) 
\dx  /,        J  !      J     0+1 

fek>\  -      (_0        +0        _0        +0        )    /    (2h3)       (2-5c) 
\dX3/.  *         J-2     ^J-l     *TJ+1         j'+2 


J 


/jl^\     =(0       -40       +6  0-40       +0       )/h4        (2-5d) 

\dx4/.  0-2      'j-l  J  j+1        j+2 

J 

in  which  h=0.1.      The   truncation  error   of  each  of  the   finite 

2 
difference   approximations    in  Equations   2-5   is    0(h    ). 

Applying  Equation  2-3   at   j=l,2,35    10,    Equations 
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2-4a  at  j=0 5  and  Equations  2 -4b  at  j=10;  and  eliminating 
the  variables  0_-,50  _,0-,-,5  and0..p  eventually  leads  to 
the  matrix  equation 
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where  I  is  the  identity  matrix.   Note  that  the  stiffness 
matrix  in  Equation  2-6  is  not  symmetric.   However^  if  the 
last  row  of  Equation  2-6  is  divided  by  two5  the  resulting 
formulation  is  symmetric.   This  results  in  the  mass  matrix 
M  no  longer  being  equal  to  I,    since  the  diagonal  element 
in  the  last  row  is  equal  to  one-half;  however^,  M  is  still 
symmetric  and  positive  definite _,  and  all  eigenvalues  are 
real.   The  symmetric  stiffness  matrix  will  be  referred  to 
hereafter  as  K  .   Note  that  all  of  the  finite  difference 
equations  used  to  construct  K  are  centered  about  j „ 


It  can  be  shown  that  this  symmetric  formulation 
is  identical  to  a  matrix  equation  developed  using  energy 
techniques  with  one -half  of  a  mass  interval  at  j=10. 
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An  alternate  matrix  formulation  of  the  governing 


difference  equations  can  be  constructed  in  a  similar  manner 


2  f 

by  using  the  0(h  )  approximation  Ref „  3 


(#).  *  ^J-3-6^.2+12^-l-10^+3^+1)/(2h3) 

(2-7) 

instead  of  Equation  2.5c. 

As  a  consequence  of  the  fact  that  Equation  2-7  is  not 
centered  about  j3  the  fictitious  node  j=12  is  not  required_, 

and  Equation  2-3  is  applied  at  j=l,253.j 9°   The 

boundary  conditions  are  satisfied  by  applying  Equations  2-4a 
at  j=0  and  Equations  2-4b  at  j=10.   Eliminating  0_,  ^  <?j>  3 
0LOj  an<^  0   leads  to  the  unsymmetric  matrix  formulation 
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(2-8) 


The  stiffness  matrix  in  Equation  2-8  will  be  referred  to 


as  K  . 
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III.   ANALYSIS  AND  RESULTS 

A  perturbation  theory  given  by  Wilkinson  Ref.  2   for 
real  symmetric  matrices  will  be  used  to  prove  that  the 
eigenvalues  of  the  real  unsymmetric  formulation  given  by 
Equation  2-8  are  real,   ^he  brief  description  of  the  theory 
which  follows  Is  expanded  in  the  Appendix. 

Let 

K   =  A  +  B 
2 

where  K  is  the  real  unsymmetric  square  matrix  of  order  n, 
A  is  an  n  x  n  real  symmetric  matrix,  and  B  is  an  n  x  n  real 
unsymmetric  matrix.   Every  eigenvalue  of  K  lies  within  at 
least  one  of  the  n  circular  discs  whose  centers  are  the 
eigenvalues  of  A.   The  radius  of  these  eirculai  discs  is 
I  B 


2,  the  2 -norm  of  B,  where 

|b||  2   =    (maximum  eigenvalue    of  B  B) 


2 


IT 

and  B  is  the  Hermitian  conjugate  matrix  of  B.   The  matrix 

H 
B  B  is  then  Hermitian,  the  eigenvalues  of  which  are  all 

real.   If  any  m  of  the  n  circular  discs  form  a  connected 

domain  isolated  from  the  others,  there  are  precisely  m 

eigenvalues  in  this  connected  domain.   In  the  case  of  real 

unsymmetric  matrices,  any  complex  eigenvalues  must  occur 

th 

as  conjugate  pairs.   If  the  i  '  circular  disc  is  isolated 

and  contains  only  one  eigenvalue,  this  eigenvalue  must 
therefore  be  real. 

The  matrix  K  can  be  decomposed  in  many  ways  to  obtain 
the  matrices  A  and  B.   One  technique  is  to  use  the  symmetric 
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matrix  K  as  A  and  obtain  B  by  subtraction.   This  method 


was  tried  and  resulted  in  such  a  large 


B  p  that  only 


two  eigenvalues  of  K  could  be  proven  real.   Consequently s 
two  other  approaches  were  considered  in  an  attempt  to 
develop  a  smaller  norm  size  that  would  prove  all  the  eigen- 
values to  be  real  and  positive.   These  two  approaches  are 
referred  to  as  Method  I  and  Method  II. 

In  Method  Ts    the  matrix  K  is  decomposed  such  that   B 
contains  the  smallest  number  of  elements  possible _,  i.e.. 
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B. 


Equation  3-1  shows  only  that  portion  of  K  that  requires 

alteration.   The  eigenvalues  of  A   and  the  2-norm  of  BT 

e  I  I 

were  determined  using  the  IBM  SSP  subroutine   EIGENP  and 
are  tabulated  in  Table  I  and  plotted  in  Figure  2.   Note 
in  Figure  2  that  five  of  the  circular  discs  are  isolated 
and  in  the  positive  domain.   As  a  result,  applying  the 
perturbation  theory  proves  that  at  least  five  of  the  nine 
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eigenvalues  are  real  and  positive.   Also  note  that  the 
five  real  eigenvalues  are  the  five  largest  eigenvalues .   In 
order  to  prove  that  the  lower  eigenvalues  in  the  overlapping 
circular  discs  are  realj,  the  inverse  formulation 


1  I  <±>  =  K"1  0 


(3-2) 


is  examined.   If  lA  can  be  proven  real  and  positive _,  then 
\  is  real  and  positive. 


TABLE  I 
EIGENVALUES  AND  NORMS  -  METHOD  I 
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AlsOj  inversion  directly  relates  the  lowest  eigenvalue  of 

Kp  to  the  largest  eigenvalue  of  K  3    and  so  on.   Since  the 

lower  eigenvalues  of  K?  have  not  yet  been  proven  real,  the 

higher  eigenvalues  of  K   are  of  particular  interest. 

The  decomposition  of  K"1  for  Method  I  is 

2 
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The  eigenvalues  of  A  (  and  the  2-norm  of  B  ,  are  also 

I  I 

shown  in  Table  I  and  in  Figure  3.   Examination  of  Figure  3 

reveals  that  the  largest  eigenvalue  of  K  "  is  real  and 

positive;  hence }    the  lowest  eigenvalue  of  Kp  is  real  and 

positive.   Thus }    six  of  the  nine  eigenvalues  have  been 

proven  real  and  positive.   To  verify  that  these  eigenvalues 

are  in  fact  real  and  positive  _,  the  actual  eigenvalues  of 

K  and  K   were  computed  using  EIGENP  and  are  tabulated 
2      2  F 

in  Table  I  and  plotted  in  Figure  2  and  Figure  3.   They  do 
in  fact  lie  within  the  circular  discs. 


--1 


;o  a 


The  decomposition  of  Kp  and  K   of  Method  I  led  t< 
relatively  large  norm  and  the  proof  that  only  six  of  the 
nine  eigenvalues  were  real  and  positive.   Consider  now 
Method  II  in  which  K   is  decomposed  such  that  the  moduli 
of  the  elements  of  B  are  as  small  as  possible  (making 
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B  skew- symmetric  )_,  I.e., 
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Only  that  portion  of  Kp  requiring  alteration  is  shown  in 

Equation  3-4.   Using  the  same  procedures  as  in  Method  I 

leads  to  the  results  tabulated  in  Table  II  and  plotted  in 

Figure  4.   The  norm  of  B   is  much  smaller  than  the  norm 

II 

of  B  .  and  examination  of  Figure  4  reveals  that  seven  of  the 
nine  eigenvalues  of  K  are  real  and  positive.   Decomposing 
K    in  a  manner  similar  to  Equation  3-4  leads  to  A   ' 
and  B   i.   The  results  are  also  shown  in  Table  II  and  in 
Figure  5.   Using  this  formulation,  the  lowest  eigenvalue 
of  Kp  is  proven  real  and  positive.   Consequently,  use  of 
Method  II  proves  eight  of  the  nine  eigenvalues  to  be  real 
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and  positive.   The  last  eigenvalue  must  also  "be  real  since 
it  can  have  no  conjugate.  It  must  also  be  positive  since  it 
is  larger  than  the  lowest  eigenvalue  which  has  been  proven 
positive.   Therefore  all  the  eigenvalues  of  K?  are  real  and 
positive . 

A  third  method  considered  is  to  treat  K  in  a  manner 
similar  to  the  conversion  of  Equation  2-6  to  the  symmetric 
formulation  K,  .   In  this  method,  K  is  altered  by  multiplying 
row  9  by  two.   Although  this  leads  to  another  unsymmetric 
matrix _,  a  much  smaller   B  ?  can  be  developed.   However, 
from  past  experience  it  is  concluded  that  this  approach 
would  not  prove  all  eigenvalues  real  and  positive  without 
resorting  to  the  inversion  of  the  matrix.    Since  Method  II 
proved  ^p  eigenvalues  of  K.  real  and  positive  without 
alteration  of  the  matrix,  and  since  the  eigenvalues  of  this 

2 

third  matrix  formulation  are  not  uu)  >    this  approach  is  not 
pursued  further. 

As  a  check  on  the  accuracy  of  the  calculated  eigenvalues 
of  Kpj  a  comparison  of  the  exact  natural  frequencies  [Ref.  4 
and  the  numerical  natural  frequencies  is  made  in  Table  III 
for  the  five  lowest  frequencies. 
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TABLE  II 
EIGENVALUES  AND  NORMS  -  METHOD  II 


i 

Eigenvalue  \: 

AH 

K2 

A         » 
AII 

v1 

1 

-0.144 

0.0013 

-5.58 

0.066 

2 

o.o4o 

0.0^9 

O.070 

0.079 

3 

0.363 

0.369 

0.084 

0.107 

4 

1.31 

1.31 

O.119 

0.168 

5 

3.17 

3.16 

0.209 

0,3l6 

6 

5.97 

5.96 

0.494 

0.764 

V 

Q ,  37 

JL  •  f  £- 

n     <7-\ 

8 

12.68 

12.68 

17.10 

20.42 

9 

15.11 

15.11 

810.6 

_ , . 

800.1 

Value  of 


II 


2 


0.148 


Value  of 


II 


2 


90.3 
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TABLE  III 
NATURAL  FREQUENCIES  OF  CANTILEVER  BEAM 


Natural  Frequenc 

y   rad/sec 

Exact 

Ref.  41 

Numerical 

O.560 

0.562 

3.510 

3.522 

9.820 

9.671 

19.20 

18.210 

31.80 

28.300 

1 
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IV.   MULTIPLE  VARIABLE  FORMULATION 

An  alternate  un symmetric  formulation  is  obtained  when 
the  two  variables  0  and  /7[  }   where  //?  is  the  bending  moment 
in  the  beam.,  are  used  instead  of  only  0  .   The  governing 
equations  are 


Ox 


and 


The  boundary  conditions  for  the  two -variable  formulation 

are 

^    A0 


[*}-  i  ) 


and 


-d^=0 


dx  at  x=1         (4-4) 

Consider  the  beam  shown  in  Figure  1.   Applying  Equation 

4-1  and  Equation  4-2  at  j=lj2,3,?.. 105  Equation  4-3 

at  j=05  Equation  4-4  at  j=10,  and  eliminating  0,0  . 
^ilJ  ^-l^^To5  and/?(7n  leads  to  the  matrix  equation 
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& 

-2 

0 

1 

0 

• 

-2 

-h2 

1 

0 

0 

•         • 

• 

0 

1 

0 

-2 

0 

1 

0    . 

• 

1 

0 

-2 

-h2 

1 

0 

0    . 

• 

0 

0 

0 

1 

0 

-2 

0 

1 

0    . 

0 

0 

1 

0 

-2 

-h2 

1 

0 

0    . 

(4-5) 
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0 

-2 
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.    0 
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0 

-2 

-h2 

1 

0 

0 
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.    0 

0 

1 

0 

-2 
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.    0 

1 

0 

-2 

-n2 

1 

•         • 

0 

0 

2 

0 

o 

0 

0 

0 
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^8 


#7 
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9 

9 


10 


^10, 


Xh' 


1 

0 

0 

0 

0 

• 

• 

0 

0 

0 

0 

0 

• 

• 
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0 

1 
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• 
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0 

0 
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0 
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• 
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• 
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0 

0 

0 

0 
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0 

0 

• 
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0 

0 
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1 

0 

0 

0 
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Equation  4-5  is  unsyrametric .   However,  if  the  variables  /77  ->} 

y?{ 2> ^flO  are  elimina"ted^  Equation  4-5  reduces 

to  Equation  2-6.   As  previously  discussed.  Equation  2-6  is 
simply  a  symmetric  matrix  in  disguise,  so  that  the  eigenvalues 
are  real  and  positive. 

The  two- variable  formulation  has  been  presented  here  to 
point  out  that  many  unsymmetric  matrices  may  in  fact  be 
reducible  to  symmetric  matrices;  hence,  they  have  all  of 
the  properties  of  symmetric  matrices . 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  basic  concept  illustrated  in  this  thesis  is  that 
for  a  given  unsymmetric  matrix  the  eigenvalues  may  be 
examined  for  proof  that  they  are  real  and  positive  without 
actually  calculating  them.   Use  of  the  perturbation  theory 
given  by  Wilkinson  gives  a  relatively  straightforward  method 
by  which  this  may  be  accomplished.   For  the  simple  example 
selected  calculation  of  the  eigenvalues  is  easily  carried  out, 
thus  confirming  the  results. 

It  is  recommended  that  this  technique  be  applied  to  other 
eigenvalue  problem  formulations  of  more  involved  physical 
systems.   If  the  technique  proves  itself  reliable  in  appli- 
cation to  matrices  of  large  order,  significant  computer  time 
can  be  saved  by  ascertaining  whether  a  given  modelling  method 
will  result  In  real,  positive  eigenvalues  prior  to  conducting 
lengthy  eigenvalue  computations. 
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APPENDIX 
PERTURBATION  THEORY  GIVEN  BY  WILKINSON 

With  slight  changes  of  notation  and.  text,  the  following 
material  is  reproduced  for  ease  of  reference  from  the  pages 
cited  for  Reference  2. 

Consider  an  eigenvalue /\  (possibly  complex)  of  the  real 

system 

(A  +  B)x  =  .Xx  (A-la) 

in  which  A  is  symmetric  but  B  is  not.   Since  A  is  real 

there  is  an  orthogonal  matrix  H  such  that 

H_1AH  -  diag  (  \.  )  (A-lb) 

and  the  X .  are  real.  Because  X  is  an  eigenvalue  of  (A-i-B), 
the  matrix  (A+B-AI)  is  singular  and  hence  its  determinant: 
is  zero.   From  Equation  A-la  and  Equation  A-lb 

H_1(A+B-Xl)H  =  diag  (\±-  \  )+H~  BH   (A-2 ) 
and  taking  determinants ,  the  matrix  on  the  right  of  Equation 
A-2  must  also  be  singular.   Two  cases  are  distinguished. 

Case  1.   \  =  \.      for  some  i 

Case  2.   X^X--   f°r  any  ^->    so  that 


.-1 


-1  -1 


diag(\i-  X)+H"  BH  =  diag(Xi-X  )  [  Miag(  \  ±-  \  )"  H  BH  j   (A-3) 

and  taking  determinants  again,  the  matrix  in  brackets  must 
be  singular.   If  (I  +Z  )9   where  Z  is  a  general  matrix,  is 


singular, 


a  !•  for  if 


J  2 


<  1,  none  of  the 
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eigenvalues  (I  +  Z  )  can  be  zero   Ref.  5 


giving 


diag  (\  -X)""^  ^H 


max  I  X.-  X|~ 
•   i    ii 


^  1 


Hence 
(A-4) 


-1 

j 

)     | 

1 

i 

H 

1 

Bl 

H 

1  2 

1     1 

2  j 

i 

2 


i.e.  (in  both  Case  1  and  Case  2) 


mm 


\-\ 

L. 

H"1 

2 

H 
i 

2 

B 

>  1   (A-5) 


(A-6) 


for  at  least  one  \.  for  every  eigenvalue  \   .   That  is  to 
say,  every  eigenvalue  of  (A  +  B)  lies  in  at  least  one  of  the 
circular  discs  given  by 


X  -  x 


< 

t   -1 

H 

p 

H 

2 

i 

B 

1 

(A-7) 


and  since  H  is  orthogonal  (all  eigenvalues  of  which  have 


f.\^^-,iTi-v-  one  '   these  vprni"^  t. o 


X.-  X 

1 


B 


(A-8) 


Further_,  if  any  s  of  these  discs  form  a  connected 
domain  isolated  from  the  others _,  then  there  are  precisely 
s  eigenvalues  in  this  domain. 

In  the  case  of  a  real  unsymmetric  matrix  such  as  (A  +  B)5 

any  complex  eigenvalues  must  occur  in  conjugate  pairs ,   If 

.th 


Xi-X, 


J 


>2 


B  2  (J  ^   i)j   then  the  i    circular  disc  is 
isolated  and  contains  only  one  eigenvalue.   This  eigenvalue 
must  therefore  be  real.   Hence,  if  all  eigenvalues  of  A 


are  separated  by  more  than  2 
(A  +  B)  are  real. 


2- 


.  all  eigenvalues  of 
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